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Abstract 

We describe a collection of computer scripts written in PARI/GP to compute, 
for reflection groups determined by finite- volume polyhedra in H'^, the commensura- 
bility invariants known as the invariant trace field and invariant quaternion algebra. 
Our scripts also allow one to determine arithmeticity of such groups and the iso- 
morphism class of the invariant quaternion algebra by analyzing its ramification. 

We present many computed examples of these invariants. This is enough to show 
that most of the groups that we consider are pairwise incommensurable. For pairs 
of groups with identical invariants, not all is lost: when both groups are arithmetic, 
having identical invariants guarantees commensurability. We discover many "unex- 
pected" commensurable pairs this way. We also present a non-arithmetic pair with 
identical invariants for which we cannot determine commensurability. 

1 Introduction 

Suppose that P is a finite- volume polyhedron in each of whose dihedral angles 
is an integer submultiple of vr. Then the group A(P) generated by reflections in 
the faces of P is a discrete subgroup of Isom(]H[^). If one restricts attention to the 
subgroup r(P) consisting of orientation-preserving elements of A(P), one naturally 
obtains a discrete subgroup of PSL{2,C) = Isom^(IH^). This very classical family 
of finite-covolume Kleinian groups is known as the family of polyhedral reflection 
groups. 

There is a complete classification of hyperbolic polyhedra with non-obtuse dihe- 
dral angles, and hence of hyperbolic reflection groups, given by Andreev's Theorem 
[3| [T6] (see also [Ml E] for alternatives to the classical proof); however, many 
more detailed questions about the resulting reflection group remain mysterious. We 
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will refer to finite-volume hyperbolic polyhedra with non-obtuse dihedral angles as 
Andreev Polyhedra and finite- volume hyperbolic polyhedra with integer submultiple 
of vr dihedral angles Coxeter Polyhedra. 

A fundamental question for general Kleinian groups is: given Ti and r2 does 
there exist an appropriate conjugating element g G PSL{2, C) so that Ti and gT2g^'^ 
both have a finite-index subgroup in common? In this case, Fi and r2 are called 
commensurable. Commensurable Kleinian groups have many properties in common, 
including coincidences in the lengths of closed geodesies (in the corresponding orb- 
ifolds) and a rational relationship (a commensurability) between their covolumes, 
if the groups are of finite-covolume. See [23] for many more interesting aspects of 
commensurability in the context of Kleinian groups. 

If Fi and r2 are fundamental groups of hyperbolic manifolds Mi and M2, com- 
mensurability is the same as the existence of a common finite-sheeted cover M of 
Ml and of M2. Similarly, if r(Pi) and r(P2) are polyhedral reflection groups, they 
are commensurable if and only if there is a larger polyhedron Q that is tiled both 
by Pi under reflections in the faces (of Pi) and by P2 under reflections in the faces 
(of P2.) The existence of such a polyhedron Q is clearly a fundamental and delicate 
question from hyperbolic geometry. See Figured] for an example of two commensu- 
rable polyhedra. These coordinates for these polyhedra were computed using [25] 
and displayed in the conformal ball model of H'^ using Geomview [1] . 

A pair of sophisticated number-theoretic invariants has been developed by Reid 
and others to distinguish between commensurability classes of general flnite-covolume 
Kleinian groups. See the recent textbook [22] and the many references therein. 
Given F, these invariants are a number fleld A;(F) known as the invariant trace field 
and a quaternion algebra j4(F) over k{T) known as the invariant quaternion algebra. 
In fact, the invariant trace fleld is obtained by intersecting all of the fields generated 
by traces of elements of the finite-index subgroups of F. It is no surprise that such a 
field is related to commensurability because the trace of a loxodromic element a of F 
is related to the translation distance d along the axis of a by 2 cosh((i) = Re(tr(j4)). 

The pair [k{r),A{T)) does a pretty good job to distinguish commensurability 
classes, but there are examples of incommensurable Kleinian groups with the same 
(A:(F), j4(F)). For arithmetic groups, however, the pair (A;(F),^(F)) is a complete 
commensurability invariant. Thus, one can find unexpected commensurable pairs 
of groups by computing these two invariants and by verifying that each group is 
arithmetic. See Subsection 17.71 for examples of such pairs that were discovered in 
this way. The precise definitions of the invariant trace field, the invariant quaternion 
algebra, and an arithmetic group will be given in Section [3l 

It can be rather difficult to compute the invariant trace field and invariant quater- 
nion algebra of a given Kleinian group "by hand." However there is a beautiful 
computer program called SNAP [10] written by Coulson, Goodman, Hodgson, and 
Neumann, as described in [9]. They have computed the invariant trace field and 
invariant quaternion algebra, as well as many other interesting invariants, for many 
of the manifolds in the Hildebrand- Weeks census [12] and in the Hodgson- Weeks 
census |15j . The basic idea used in SNAP is to compute a high-precision decimal 
approximation for an ideal triangulation of the desired manifold M using Newton's 



2 




Figure 1: Two commensurable polyhedra Pi (left top) and P2 (right top) which tile a 
common larger polyhedron Q, here the right-angled dodecahedron. 
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Method and then to use the LLL algorithm [19] to guess exact algebraic numbers 
from the approximate values. These guessed values can be checked for correctness 
using the gluing equations describing M, and if the values are correct, the invari- 
ant trace field and invariant quaternion algebra can be computed from the exact 
triangulation. 

SNAP provides a vast source of examples, also seen in the appendix of the book 
[22], and adds enormous flavor to the field. The fundamental techniques used in 
SNAP provide inspiration for our current work with polyhedral reflection groups. 

In the case of polyhedral reflection groups there is a simplified description of 
the invariant trace field and the invariant quaternion algebra in terms of the Gram 
matrix of the polyhedron [21]. This theorem avoids the rather tedious trace calcu- 
lations and manipulation of explicit generators of the group. Following the general 
technique used in the program SNAP we compute a set of outward unit normals 
to the faces of the polyhedron P to a high decimal precision and then use the LLL 
algorithm to guess the exact normals as algebraic numbers. From these normals the 
Gram matrix is readily computed, both allowing us to check whether the guessed 
algebraic numbers are in fact correct, and providing the exact data needed to use 
the theorem from [21] in order to compute the invariant trace field and invariant 
quaternion algebra for T(P). 

Our technique is illustrated for a simple example in Section [5] and a description 
of our program (available to download, see [4]) is given in Section [6l Section |4] 
provides details on how to interpret the quaternion algebra. Finally in Section [7] we 
provide many results of our computations. 

Acknowledgments: We thank Colin Maclachlan and Alan Reid for their beautiful 
work and exposition on the subject in [22]. We also thank Alan Reid for his many 
helpful comments. 

We thank the authors of SNAP [10] and the corresponding paper [9], which 
inspired this project (including the choice of name for our collection of scripts). 
Among them Craig Hodgson has provided helpful comments. 

We thank Andrei Vesnin who informed us about his result about arithmeticity 
of Lobell polyhedra. 

We effusively thank the writers of PARI/GP [23, the system in which we have 
written our entire program and which is also used in SNAP |10j . 

The third author thanks Mikhail Lyubich and Ilia Binder for their financial sup- 
port and interest in the project. He also thanks John Hubbard, to whom this volume 
is dedicated, for introducing him to hyperbolic geometry and for his enthusiasm for 
mathematics in general and experimental mathematics in particular. 

2 Hyperbolic polyhedra and the Gram Matrix 

We briefly recall some fundamental hyperbolic geometry, including the definition of 
a hyperbolic polyhedron and of the Gram matrix of a polyhedron. 
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Let E^'^ be the four-dimensional Euclidean space with the indefinite metric 
||x|p = —Xq + Xi + X2 + X3. Then hyperbolic space is the component having 
xq > of the subset of E^'^ given by 

||x|| — Xq ~|~ X-^ ~\~ ~\~ 37^ — X 

with the Riemannian metric induced by the indefinite metric 

— + dx\ + dx^ + dx\. 

The hyperplane orthogonal to a vector v G E'^'^ intersects 'M? if and only if 
(v, v) > 0. Let V G E^'^ be a vector with (v, v) > 0, and define 

= {w G El3|(w,v) = 0} 
to be the hyperbolic plane orthogonal to v; and the corresponding closed half space: 

H+ = {^e H3|(w,v) > 0}. 

It is a well known fact that given two planes and P-^ in with (v, v) = 1 and 
(w,w) = 1, they: 

• intersect in a line if and only if (v, w)^ < 1, in which case their dihedral angle 
is arccos(— (v, w)). 

• intersect in a single point at infinity if and only if (v,w)^ = 1, in this case 
their dihedral angle is 0. 

• are disjoint if and only if (v, w)^ > 1, in which case the distance between them 
is arccosh(— (v, w))). 

Suppose that ei, . . . ,e„ satisfy (ei,ej) > for each i. Then, a hyperbolic poly- 
hedron is an intersection 

n 
1=0 

having non-empty interior. 

If we normalize the vectors that are orthogonal to the faces of a polyhedron 
P, the Gram Matrix oi P is given by Mij{P) = 2(ej,ej). It is also common to define 
the Gram matrix without this factor of 2, but our definition is more convenient for 
arithmetic reasons. By construction, a Gram matrix is symmetric and has 2s on 
the diagonal. Notice that the Gram matrix encodes information about both the 
dihedral angles between adjacent faces of P and the hyperbolic distances between 
non-adjacent faces. 
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3 Invariant Trace Field, Invariant Quaternion 
Algebra, and Arithmeticity 



The trace field of a subgroup T of PSL{2, C) is the field generated by the traces of its 
elements; that is, Q(trr) := Q(tr7 : 7 G r)0 This field is not a commensurablitity 
invariant as shown by the following example found in |22j . 
Consider the group T generated by 



It is easy to see that X normalizes F and its square is the identity (in PSL{2, C)), so 
that A = (r, X) contains T as a subgroup of index 2 and is therefore commensurable 



contains i in addition to uj. 

The easiest way to fix this, that is, to get a commensurability invariant related 
to the trace field, is to associate to T the intersection of the trace fields of all finite 
index subgroups of F; this is the invariant trace field denoted kT. 

While this definition clearly shows commensurability invariance, it does not lend 
itself to practical calculation. The proof of Theorem 3.3.4 in [22] brings us closer: 
it shows that instead of intersecting many trace fields, one can look at a single 
one, namely, the invariant trace field of F equals the trace field of its subgroup 
r(2) — . J ^ Y). (This also shows that the invariant trace field is non-trivial 
which is not entirely clear from the definition as an intersection.) When a finite set of 
generators for the group is known, this is actually enough to compute the invariant 
trace field. Indeed, Lemma 3.5.3 in [22j establishes that if F = (71,72, • • • ,'Jn), the 
invariant trace field of F is generated by {tr(7j) : 1 < i < n} U {tr(7j7j) : 1 < i < 
j <n}U {tr(7i7j7fc : I < i < j < k < n}. 

For reflection groups a more efficient description can be given in terms of the 
Gram matrix. The description given above uses roughly generators for a 

polyhedron with n faces; the following description will use only around v? 12. But 
before we state it we need to define a certain quadratic space over the field k{P^ := 
: {n, «2, • • • , v} C {1, 2, . . . , n}) associated to a polyhedron; this 
space will also appear in the next section in the theorem used to calculate the 
invariant quaternion algebra. 

As in the previous section, given a polyhedron P we will denote the outward- 
pointing normals to the faces by ei , . . . , e„ and the Gram matrix by (ojj ) . Define 
M{P) as the vector space over k{P) spanned by of all the vectors of the form 
Gi^ where {ii, 12, . . . , v} ranges over the subsets of {1,2,..., n} 
and n is the number of faces of P. This space M{P) will be equipped with the 
restriction of the quadratic form with signature (3, 1) used in 'M? . We recall that 

^Note that for 7 G PSL{2,C), the trace tr7 is only defined up to sign. 




where w = (-l + iV3)/2. The trace field of F is Q{y/^). Now let X 





trace field of A 
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the discriminant of a non-degenerate symmetric bilinear form (•, •) is defined as 
det i{vi,Vj))-j where {vi}i is a basis for the vector space on which the form is defined. 
The discriminant does depend on the choice of basis, but for different bases the 
discriminants differ by multiphcation by a square in the ground field: indeed, if Ui = 
OiijVj, the discriminant for the basis {ui}i is that of the basis {vi}i multiplied 
by det(Qij)2. 

Now we can state the theorem we use to calculate the invariant trace field, 
Theorem 10.4.1 in [22] : 

Theorem 1 Let P be a Coxeter polyhedron and let T be the reflection group it 
determines. Let (oij) be the Gram matrix of P. The invariant trace field ofT is 
k{P)[\fd) , where d is the discriminant of the quadratic space M{P) and k{P) = 
Qioij^i^ai^i^ ■ ■ ■ ai^i^ : {ii,i2, ■ ■ ■ ,ir} C {1,2, . . . is the field defined previously. 

3.1 The Invariant Quaternion Algebra 

A quaternion algebra over a field F is a four-dimensional associative algebra A with 
basis {l,i,j, /c} satisfying = al, f' = bl and ij = ji = —k for some a,b ^ F. 
Note that A;^ = (ij)^ = ijij = —ijji = —abl. The case F = M, a = 6 = — 1 gives 
Hamilton's quaternions. 

The quaternion algebra defined by a pair a,b of elements of F is denoted by its 
Hilbert symbol A quaternion algebra does not uniquely determine a Hilbert 

symbol for it, since, for example, = ^ "'""^ ^ = ^ fg^ any invertible 

elements u,v £ F. Fortunately, there is a computationally effective way of deciding 
whether two Hilbert symbols give the same quaternion algebra. This will be dis- 
cussed in Sectional for now we will just define the invariant quaternion algebra of a 
subgroup of PSL{2, C) and state the theorem we use to calculate a Hilbert symbol 
for it. 

Given any non-elementar}il subgroup F of PSL{2,C), we can form the algebra 
:= {"^CLi'ji : Qi G Q(trF),7j G F}. (Abusing notation slightly we consider the 
elements of F as matrices defined up to sign.) This turns out to be a quaternion 
algebra over the trace field Q(trF) (see Theorem 3.2.1 in |22j). 

Just as with the trace fields, we define the invariant quaternion algebra of F, 
denoted by ^F, as the intersection of all the quaternion algebras associated to 
finite-index subgroups of F. 

When F is finitely generated in addition to non-elementary, we are in a situation 
similar to that of the invariant trace field in that the invariant quaternion algebra is 
simply the quaternion algebra associated to the subgroup F(2) of F, or in symbols, 
AT = AoF(2). To see this, note that Theorem 3.3.5 in [22] states that for finitely 
generated non-elementary F, the quaternion algebra ^or^^"* is a commensur ability 
invariant. Now, given an arbitrary finite- index subgroup A of F we have ^F C 
AoF(2) = AoA(2) C AqK. 

^This means that the action of F on H'^ U C has no finite orbits. Reflection groups determined by 
finite- volume polyhedra are always non-elementary. 
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Figure 2: The possible vertices in the singular locus of a 3-dimensional orbifold. The last 
three have vertex stabilizer containing /I4. 

In the case where T is the reflection group of a polyhedron P, AT can be identified 
as the even-degree subalgebra of a certain Clifford algebra associated with P. Let us 
briefly recall the basic notions related to Clifford algebras. Given an n-dimensional 
vector space V over a fleld F equipped with a non-degenerate symmetric bilinear 
form (•, •) and associated quadratic form ||-||^, the Clifford algebra it determines 
is the 2^-dimensional F-algebra generated by all formal products of vectors in V 
subject to the condition v'^ = {v,v) 1 (where 1 is the empty product of vectors). If 
{vi,V2, . . . , Vn} is an orthogonal basis of V, a basis for the Clifford algebra C{V) is 
{viiVi^ • ■ ■ Vi^ : 1 < r < n,l < ii < 12 < • • • < ir < n}. There is a Z2-grading of 
C{V) given on monomials by the parity of the number of vector factors. 

Now we can state the result mentioned above: AT is the even-degree subalgebra 
of C (M) where M is the vector space over the invariant trace fleld of T that appears 
in Theorem [TJ This relationship between the invariant trace algebra and M allows 
one to prove a theorem giving an algorithm for computing the Hilbert symbol for 
the invariant quaternion algebra, part of Theorem 3.1 in [21] : 

Theorem 2 The invariant quaternion algebra of the reflection group T of a poly- 
hedron P is given by 




where {ui,U2,U3,U4} is an orthogonal basis for the quadratic space M{P) defined 
in the previous section. 

In many cases when studying an orbifold O = M^/T, a simple observation about 
the singular locus of O leads to the fact that the invariant quaternion algebra can 
represented by the Hilbert symbol ^ "^.'^T^ ^ . This happens particularly often for 
polyhedral reflection groups. 

Any vertex in the singular locus must be trivalent and must have labels from 
the short list shown in Figure [2l See [5] for more information on orbifolds, and in 
particular page 24 from which Figure [2] is essentially copied. 

The last three have vertex stabilizer containing ^4, so if H'^/F has singular 
locus containing such a vertex, F must contain as a subgroup. In this case, the 
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invariant quaternion algebra can be represented by the Hilbert symbol |^ ^'p j , see 

[22], Lemma 5.4.1. (See also Lemma 5.4.2.) 

For a polyhedral reflection group generated by a Coxeter polyhedron P, the 
corresponding orbifold H^/r has underlying space and the singular set is (an 
unknotted) copy of the edge graph of P. The label at each edge having dihedral 
angle ^ is merely n. In many cases, this simplification makes it easier to compute 
the Hilbert symbol of a polyhedral reflection group. We do not automate this check 
within our program, but it can be useful to the reader. 



4 Quaternion algebras and their invariants 

As mentioned before, a quaternion algebra is fully determined by its Hilbert symbol 
, although this is by no means unique. For example, 

b,a\ fa,—ab\ ( ax^ ^by^ 

and 



F J \ F J ' \ F 

all determine the same algebra (here x and y are arbitrary invertible elements of 

Taking F = M, the Hilbert symbol ( "'''m""^ ) represents the ordinary quaternions 



(or Hamiltonians,) denoted by Ti. If F is any field, the Hilbert symbol y^j is 

isomorphic to M2{F), the two-by-two matrices over F. 

A natural question now arises; namely, when do two Hilbert symbols represent 
the same quaternion algebra? This question is pertinent for us — especially in the 
case when F is a number field — because the invariant quaternion algebra is an 
invariant for a refiection group. For these purposes, we will need some way of 
classifying quaternion algebras over number fields. All of the following material on 
quaternion algebras appears in the reference |22j . 

The first step in the classification of quaternion algebras is the following theorem: 



Theorem 3 Let A be a quaternion algebra over a field F . Then either A is a 
division algebra or A is isomorphic to M2{F). 

In the latter case, we say that A splits. There are several different ways of expressing 
this condition, one of which will be particularly useful for us: 

Theorem 4 The quaternion algebra A = splits over F if and only if the 

equation 

ax^ + = 1 (1) 
has a solution in F x F. We call this equation the Hilbert equation of A. 
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When F is a number field, it turns out that in order to classify the quaternion 
algebras over F completely we need to look at quaternion algebras over the com- 
pletions of F with respect to its valuations. The first chapter of [22] contains a 
brief introduction to number fields and valuations, and [18] is a standard text on 
the subject. 

Definition 5 Let F he any field. A valuation on F is a map v : F ^ such that 
(i) u{x) > for all X £ F and z^(x) = if and only if x = 0, 
(a) u{xy) = u{x)v[y) for all x,y £ F, and 
(Hi) v{x + y) < v{x) + iy{y) for all x,y £ F. 

Any field admits a trivial valuation z/(x) = 1 for all x ^ 0. When is a subfield 
of the real (or complex) numbers, the ordinary absolute value (or modulus) function 
is a valuation when restricted to F. In general, valuations on a field fall into two 
different classes. 

Definition 6 If a valuation v on a field F also satisfies 
(iv) ^{x + y) < max{i^{x), I'iy)} for all x,y £ F, 

then V is called a non-Archimedean valuation. If the valuation v does not satisfy 
(iv), then it is called Archimedean. 

There is also a notion of equivalence between valuations. 

Definition 7 Two valuations v\ and V2 on F are called equivalent if there exists 
some a £ M"*" such that V2{x) = {i'i{x))" for all x £ F. 

When F is a number field, it is possible to classify all valuations on F up to 
this notion of equivalence. Let cr be a real or complex embedding of F. Then a 
valuation fo- can be defined by fo-(x) = |ct(x)|, where | • | is the absolute value on M 
or modulus on C. This is an Archimedean valuation on F, and up to equivalence 
these are the only Archimedean valuations that F admits. 

Denote by Rp the ring of integers of F — i.e. the set of elements of F satisfying 
some monic polynomial equation with integer coefficients — which is a subring of F. 
Let V he a, prime ideal in Rp. Define a function n-p : Rp ^ Z hy n-p{a) = m, 
where m is the largest integer such that a £ V"^. Since F is the field of fractions 
of Rp, n-p can be extended to all of F by the rule n-p{a/b) = np{a) — np{b). Now 
pick c with < c < 1. The function u-p : Rp ^ M"*" given by i>p{x) = c^'p^^'^ is a 
non- Archimedean valuation on F. Moreover, all non- Archimedean valuations on F 
are equivalent to a valuation of this form. 

We summarize these facts in the following theorem. 

Theorem 8 Let F he a number field. Then every Archimedean valuation of F 
is equivalent to Vfj for some real or complex embedding a of F, and every non- 
Archimedean valuation of F is equivalent to up for some prime ideal V of Rp. The 
former are sometimes called infinite places, while the latter are called finite places. 
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A valuation u on a field F defines a metric on F by d{x,y) = u{x — y). The 
completion of F with respect to this metric is denoted by F^. Equivalent valuations 
give rise to the same completions. \{ u = for a real or complex embedding a of 
F, then Fy is isomorphic to M or C respectively. 11 v = v-p for some prime ideal 
V Q Rp, then F^, is called a V-adic field. 

Let A be a quaternion algebra over a number field F, and let F^, be the com- 
pletion of F with respect to some valuation u. Then we can construct the tensor 
product A (dp F^, which turns out to be a quaternion algebra over F^. Indeed, 
if ^ = then A (S>f F^, = Ultimately, the classification of quaternion 

algebras over F will be reduced to the classification of quaternion algebras over 
completions of F with respect to valuations u. The two following theorems will be 
useful in this regard. 

Theorem 9 Let M be the real number field. Then Ti. is the unique quaternion 
division algebra over R. 

Theorem 10 Let Fy be a V-adic field. Then there is a unique quaternion division 
algebra over Fy. 

Thus when ^4 is a quaternion algebra over a number field F and Ay is the 
corresponding quaternion algebra over a real or P-adic completion Fy of F, by 
Theorem [3] there are two possibilities: A is the unique quaternion division algebra 
over Fy, or A = M2{Fy). In the former case, we say that A ramifies at z/, while in 
the latter case we say that A splits at v. When Fy = M, there is a simple test to 
determine which algebra is represented by the Hilbert symbol ( ) : 



Various tests exist for Vy , but often the simplest test is to determine if there exists 
a solution to equation ([T]). (See Appendix El) Notice that by Theorem H] every 
quaternion algebra over the complex numbers is isomorphic to M2(C). 

The following theorem provides the necessary criterion for distinguishing be- 
tween quaternion algebras over number fields. 

Theorem 11 (Vigneras J 31^.) Let F be a number field. For each quaternion al- 
gebra A over F , denote by Ram[A) the set of all real or finite places at which A 
ramifies. Then two quaternion algebras A and A' over F are equal if and only if 
Ram{A) = Ram{A'). 

Thus the complete identification of a quaternion algebra A over a number field 
F amounts to determining Ram{A). It is easy to check if A ramifies at the real 
infinite places of F. Let a be a primitive element of F (i.e. an element whose powers 
form a basis for F over Q.) Every embedding of -F in M can be obtained from a real 
root Ui of the minimal polynomial of a over Q by extending the map fij : q — > 
linearly to F. If ^4 = and a and b are expressed as polynomials in a, then it 

is straightforward to check if condition ([2]) holds for (Ti{a) and ai{b). 




if a and b are both negative, then A = 7i, otherwise A splits. 



(2) 
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The finite places of F are more difficult to check. The most straightforward 
method is to check to see if the Hilbert equation ([T|) has a solution; our procedure 
for doing this is detailed in Appendix Rl 

4.1 Arithmeticity 

The notion of an arithmetic group comes from the theory of algebraic groups and is 
a standard way of producing finite-covolume discrete subgroups of semi-simple Lie 
groups. To see how this general theory relates to Kleinian groups, see [13j or |22j . 

In the case of Kleinian groups, the following definition coincides with the most 
general one, and is naturally related to the quaternion algebras which we have 
already mentioned. 

Let ^ be a quaternion algebra over a number field F and denote hy Rp the ring 
of integers in F. An order O in ^4 is an i^^-lattice (spanning A over F) that is 
also a ring with unity. For every complex place v F there is an embedding of 
A — > M2(C) determined by the isomorphism A®p F^ = M2(C). Given a complex 
place and an order O, we can construct a subgroup of SL{2,C), and hence of 
P5L(2,C), by taking the image Yq of the elements of O with unit norm under the 
embedding A — > M2(C) defined above. 

In the case that F has a unique complex place z> and that A ramifies over every 
real place of F, then Tq ■= Tq is a discrete subgroup of PSL{2,C) (see Sections 
8.1 and 8.2 of [22]). 

Definition 12 A Kleinian group T is called arithmetic if it is commensurable with 
To for some order O of a quaternion algebra that ramifies over every real place and 
is defined over a field with a unique complex place. 

Viewing Isom^(E['^) as S'0"'"(3, 1) furnishes an alternative construction of arith- 
metic Kleinian groups as follows. Let F be a real number field, and {V, q) a 
four-dimensional quadratic space over F with signature (3,1). Any F-linear map 
a : V — > V preserving q can be identified with an element of SO~^{3, 1), and thus 
of Isom'^(EI'^), by extension of scalars from F to M. Given an RF-laitice L C V of 
rank 4 over Rp, the group SO{L) := {a G 50+(3,l) n GL{A,F) : a{L) = L} is 
always discrete and arithmetic. 

Moreover, the groups of the form SO{L) give representatives for the commen- 
sur ability classes of all Kleinian groups that possess a non-elementary Fuchsian 
subgroup. Since all reflection groups determined by finite-volume polyhedra have 
non-elementary Fuchsian subgroups, for our purposes this can be considered the 
definition of arithmeticity. See [131 page 143] for a discussion. The distinction be- 
tween arithmetic groups arising from quaternion algebras and those arising from 
quadratic forms is also discussed in [331 pages 217-221], whose authors call the 
latter "arithmetic groups of the simplest kind". 

Aside from its relationship to algebraic groups, arithmeticity is interesting for 
many reasons including the fact that for arithmetic groups F the pair (kV, AT) is 
a complete commensurability invariant (see Section 8.4 of [22]). This will allow us 
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to identify several unexpected pairs of commensurable reflection groups which are 
presented in Section [TTTl 

To decide whether a given reflection group T determined by a polyhedron P is 
arithmetic there is a classical theorem due to Vinberg |32j : 

Theorem 13 Let {aij) he the Gram matrix of a Coxeter polyhedron P. Then the 
reflection group determined by P is arithmetic if and only if the following three 
conditions hold: 

1. K := Q{aij) is totally real. 

2. For every embedding a : K — > C such that cr|A;(p) 7^ id (where k{P) is the 
field defined in Theorem\^, the matrix {a{aij)) is positive semi-definite. 

3. The aij are algebraic integers. 

More generally, for any finite-covolume Kleinian group F, Maclachlan and Reid 
have proved a similar result, Theorem 8.3.2 in [22J: 

Theorem 14 A finite-covolume Kleinian group T is arithmetic if and only if the 
following three conditions hold: 

1. kT has exactly one complex place. 

2. AT ramifies at every real place of kT. 

3. trj is an algebraic integer for each 7 e T. 

5 Worked example 

A Lambert cube is a compact polyhedron realizing the combinatorial type of a cube, 
with three disjoint non-coplanar edges chosen and assigned dihedral angles j, ^, 
and ^, and the remaining edges assigned dihedral angles ^. It is easy to verify that 
if l,m,n > 2, then, such an assignment of dihedral angles satisfies the hypotheses of 
Andreev's Theorem. The resulting polyhedron is called the (/, m, n)-Lambert Cube, 
which we will denote by Pi^rn,n- 

In this section we illustrate our techniques by computing the invariant trace field 
and invariant quaternion algebra associated to the (3, 3, 6) Lambert cube. 

The starting point of our computation is a set of low-precision decimal approxi- 
mations of outward-pointing normal vectors {ei, . . . , eg} to the six faces of our cube. 
For a given compact hyperbolic polyhedron, it is nontrivial to construct such a set 
of outward-pointing normal vectors. One way is to use the collection of Matlab 
scripts described in [25]. Throughout this paper we will always assume the follow- 
ing normalization for the location of our polyhedron: the first three faces meet at a 
vertex, the first face has normal vector (0, 0, 0, *), the second face has normal vector 
(0, 0, *, *), and the third has form (0, *, *, *), where * indicates that no condition is 
placed on that number. 

We then use Newton's Method with extended-precision decimals to improve this 
set of approximate normals until they are very precise. (Here, we do this with 
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precision 40 numbers, but we display fewer digits for the reader.) The vectors 
{ei, . . . eg} are displayed as rows in the following matrix: 



0.0 
0.0 

-0.0 



-0.0 
0.0 

-1.00000002377892 



1.38941010745090 0.86602538319131 
0.79708547435960 1.27880618290479 
0.62728529885922 0.0 



0.0 -0.99999996237702 

0.86602540463740 0.50000001881149 

-0.0 1.0 X lO"^'' 

-0.73831913868376 1.27880621354777 

0.0 -0.0 

-1.18046043888844 1.0 x 10"^^ 



The normalization we have chosen for the location of our polyhedron assures us 
that each of these decimals should approximate an algebraic number of some (low) 
degree. There are commands in many computer algebra packages for guessing the 
minimal polynomial that is most likely satisfied by a given decimal approximate. 
Most of these commands are ultimately based on the LLL algorithm [19j . (We have 
used the command minpolyO in Maple, the command algdepO in Pari/GP, and 
the command RootApproximant [] in Mathematica 6.) Each of these commands 
requires a parameter specifying up to what degree of polynomials to search. In this 
case we specify degree 30. The resulting matrix of guessed minimal polynomials is: 



X X 
X X 
X l + X 

-9-88X2 + 48X^ -3 + 4X2 

-9 + 4X2 + 16X4 3-28X2 + 16X4 



-1-X2 + 9X^ 



X 



X 

-3 + 4X2 
X 



l + X 
-1 + 2X 
X 



1 - 28X2 + 48X4 3-28X2 + 16X4 
X X 
-3 + X + 3X2 X 



The next step is to specify which root of each given minimal polynomial is 
closest to the decimal approximate above. In the current case, the solutions of 
each polynomial are easily expressed by radicals, so we merely pick the appropriate 
expression. (For more complicated examples, our computer program uses a more 
sophisticated way of expressing algebraic numbers as described in Section [H) For 
the current example, the following matrix contains as rows our guessed exact values 
for {ei, . . .eg}. 



14 
















-1 












2 


1 
2 


N : = 







-1 








^33+6 1 


/37 




-V42+6\/37 


Vl4+2x/37 




6 




2 


12 


4 




V^-2+2 
4 




a/i4+2 V37 
4 










\/2+2v 
L 6 







6 






Corresponding to this set of guessed normal vectors we have the Gram Matrix 
Gij = 2{ei,ej): 



G := 



2 


-1 





-y/u+2V37 
2 








-1 


2 











-v^(l+\/37) 
6 








2 




-V14+2V37 

2 





\/l4+2 v/37 
2 







2 














-^14+2^37 





2 


-1 


2 





-v^(l+V37) 








-1 


2 



By checking that there are 2's down the diagonal of G and that there is —2 cos(Qij) 
in the ij-th. entry of G if faces i and j are adjacent, we can see that the guessed 
matrix N whose rows represent outward-pointing normal vectors was correct. Here 
there are 4 • 6 — 6 = 18 equations that we have checked, consistent with the number 
of guessed values in the matrix N. Consequently G is the exact Gram matrix for 
the (3, 3, 6) Lambert cube. We can use this to compute the invariant trace field and 
the invariant quaternion algebra. 

The nontrivial cyclic products from G correspond to non-trivial cycles in the 
Coxeter symbol for ^3,3,6) which is depicted in Figure[3]with the appropriate element 
of the Gram matrix written next to each edge. (For those unfamiliar with Coxeter 
symbols, seeEl) 

The non-trivial cyclic products correspond to closed loops in the Coxeter symbol. 
Always included are the squares of each entry of G, of which the only two irrational 
ones are: ((712)^ = '^^'^ and (535)^ = ^^^^^ . The other nontrivial cyclic product 
corresponds to the closed loop in the Coxeter symbol: gi25265655'535'345'4i = 11 + 
2^/37. 

Thus, k{P3,3fi) = Q(\/37). 

Notice that vi := ei, V12 := 51262 = -62, V143 := 51454363 = ^^^"^^^^ • \/3e3, 
and vi4 := 51464 = ~Vi4+2\/ 37^^ linearly independent, so that they span the 
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.956 = -1 

P o 




Figure 3: Coxeter symbol for the (3,3,6) Lambert cube. Next to each edge we display 
the corresponding element of the Gram matrix G. 



quadratic space (M, q). (See the description of Theorems [J and [2l) We now compute 
the matrix representing q with respect to the basis {vi, vi2, V143, v^}. 



7+V37 



2 10 
12 

21 +3\/37 21+3y^ 



7+^ ^i±|^ 7 + V37 



The determinant of this matrix, and hence the discriminant of the quadratic form 
q (a number well-defined up to a square in the field A;(P3 3 6)) is d = 2973-489v^ ^ 



consequently kTs^3^Q = A:(P3,3,6 
we actually have A;r3^3^6 = 



2973-489^37 \ _ g-^^^g ^ -g primitive for kiP^^s^e), 
2973-489V37 \ ^ rpj^-^ gxpressiou stih looks rather 



cumbersome, and by writing a minimal polynomial for 



2973-489V37 



and using the 



"polredabsO" command in Pari, we can check that \/—10 — 2 \/37 also generates 
this field, hence A;r3 3.6 = Q (V-W-2y/37] 



In order to use Theorems [T] and [2] to compute the invariant quaternion alge- 
bra we need to express the quadratic form q with respect to an orthogonal basis 
{wi, v^2 5 W3, W4}. The result is: 
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2 

I 



21 + 3V37 









'151-25^37 
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Thus, a Hilbert symbol describing AT^^^^q is given by 



-g(wi)9(w2),-g(wi)<;(w3) \ 
^^3,3,6 ) ' 



-3,-42-6V37 



That is, ^r3,3,6 = 1 / , 1 . 

\ Q(V-10-2y37j / 

Without using the machinery described in Section S] is difficult to interpret this 
Hilbert symbol. The invariant trace field fera^s^g has two real places, and ^Ta^s g is 
ramified at each of these two places. By attempting to solve the Hilbert Equation 
(see Appendix we also observe that 3^5 is ramified over exactly two finite 
prime ideals in the ring of integers from /cFs^a^e- These prime ideals lie over the 
rational prime 3 and we denote them by and V'^- Thus, according to Vigneras 
[31j . this finite collection of ramification data provides a complete invariant for the 



isomorphism class of AT 



3,3,6 




-10-2 V37j 

Using either Theorem 1131 or Theorem 1141 one can also deduce that Ta s^g is not 

V14+2. 



arithmetic because the Gram matrix contains the element 



37 



which is not 



an algebraic integer. However the other two hypotheses of each of these theorems 
are satisfied. 



6 Description of the program 

We've written a collection of PARI/GP scripts that automate the procedure from 
Section [5] for (finite-volume) hyperbolic polyhedra. The scripts take as input two 
matrices: the matrix of face normals, whose rows are low-precision decimal approxi- 
mations to the outward-pointing normal vectors (normalized in as in the previous 
section); and the matrix of edge labels, a square matrix ijiij) whose diagonal entries 
must be one, and whose off-diagonal terms n^j describe the relation between the i-th 
and j-th faces: n^j = means the faces are non-adjacent and any other value means 
that they meet at a dihedral angle of Txjruj. The matrix of edge labels actually 
determines the polyhedron uniquely (up to hyperbolic isometry); one can use 
to compute the approximate face normals from the edge labels. 

In a typical session, after loading the input matrices, the user runs Newton's 
Method to obtain a higher-precision approximation for the face normals, and then 
can have the computer guess and verify exact values for the normal vectors and for 
the Gram matrix. With the exact Gram matrix, the user can request the invariant 
trace field (described by a primitive element) and the invariant quaternion algebra 
(described by a Hilbert symbol, or, after an additional command, by ramification 
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data). There is also a function to test the polyhedron for arithmeticity — this too 
takes the Gram matrix as input. 

The scripts are available at ; the package includes a sample session and a user 
guide. The scripts, as mentioned before, are written in the high-level language GP, 
which helped in our effort to make the source code as readable as possible. In fact, 
the reader can consider the source code as executable statements of the theorems 
in the previous sections. 

There are two main technical challenges involved in writing these scripts: choos- 
ing an appropriate representation of general algebraic numbers, and choosing a 
systematic way of listing the non-zero cyclic products ai^i^ai^i-^ ■ ■ ■ ai^i^ from the 
Gram matrix. 

We describe a given algebraic number a as a pair (p(z),q), where p{z) is the 
minimal polynomial for a over Q and d is a decimal approximation to a. This 
representation is rather common — it is described in the textbook on computational 
number theory by Cohen [7]. While this representation is not already available in 
PARI/GP [27], it is easy to program an algebraic number package working within 
PARI/GP for this representation. 

A disadvantage of this representation for algebraic numbers is that when per- 
forming arithmetic on algebraic numbers, it is usually necessary to find composita 
of the fields generated by each number. This is not only a programming difficulty, 
it is the slowest part of our program. In some cases, when we know that we will do 
arithmetic with a given set of numbers {qi , . . . , Onji we are able to speed up our cal- 
culations by computing a primitive element (3 for Q (ai, . . . ,Q;n) and re-expressing 
each Oj as an element of Q[z\/ {p{z)) where p{z) is the minimal polynomial for jS 
over Q. (PARI/GP has a data type called "polmod" for this representation.) Once 
expressed in terms of a common field, algebraic computations in terms of these 
polmods are extremely fast in PARI/GP. 

Of course we don't typically have a priori knowledge of the outward-pointing 
normal vectors or of the entries in the Gram matrix for our polyhedron P. However, 
decimal approximations can be obtained using [25]. Just like in Section [U we use 
the LLL algorithm [19] to guess minimal polynomials for the algebraic numbers 
represented by these decimal approximations. Typically a rather high-precision 
approximation is needed (sometimes 100 digits of precision) to obtain a correct 
guess. In this case, we use Newton's Method and the high-precision capabilities 
of PARI/GP to improve the precision of the normal vectors obtained from |25| . 
Correctness of the guess is verified once we use the guessed algebraic numbers for 
the outward-pointing normals to compute the Gram matrix and verify that each 
entry corresponding to a dihedral angle ^ has the correct minimal polynomial for 
—I cos (^). Additionally we check that the diagonal entries are exactly 2. In some 
cases the guessed polynomials are not correct, but typically, by sufficiently increasing 
the number of digits of precision, reapplying Newton's Method, and guessing again, 
we arrive at correct guesses. 

Note that the number of checks, i.e. equations, equals the number of variables. If 
P has n faces, there are 4n — 6 variables: one for each coordinate of each outward- 
pointing normal, minus the six coordinates normalized to 0. Since we deal with 
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compact Coxeter polyhedra and these necessarily have three faces meeting at each 
vertex, the number of edges is 3n — 6; there is one equation for each of these, and 
there are n additional equations for the diagonal entries of the Gram matrix, giving 
a total of 4n — 6 equations. 

The guess and check philosophy is inspired by SNAP |T0], which also goes 
through the process of improving an initial approximation to the hyperbolic struc- 
ture, guessing minimal polynomials, and verifying. Such techniques are also central 
to many areas of experimental mathematics in which the LLL algorithm is used to 
guess linear dependencies that are verified a posteriori. 

When it comes to listing the non-zero cyclic products from the Gram matrix, 
we try to avoid finding more of them than necessary to compute the field they 
generate. It is easier to discuss these products in terms of the Coxeter symbol for 
the polyhedron. Recall that the Coxeter symbol is a graph whose vertices are the 
faces of the polyhedron with edges between pairs of non-adjacent faces and also 
between pairs of adjacent non-perpedicular faces. Typically edges between non- 
adjacent faces are drawn dashed, and edges between adjacent ones meeting at an 
angle of vr/n are labeled n — 2. A non-zero cyclic product in the Gram matrix 
corresponds to a closed path in the Coxeter symbol. 

An example Coxeter symbol appears in Figure [3l While this Coxeter symbol and 
most that appear in the literature are planar, this is not typically the case due to the 
large number of dashed edges. The reason so many planar Coxeter symbols appear 
in the literature is that they are especially useful for tetrahedra in high dimensions 
and those are planar. 

Our method, then, is to list all the squares of the elements of the Gram matrix, 
one for each edge in the Coxeter symbol, and a set of cycles that form a basis for 
the Z/2Z-homology of the Coxeter symbol. These allow one to express any cyclic 
product in the Gram matrix as a product of a number of basic cycles and squares or 
inverses of squares of Gram matrix elements corresponding to edges traversed more 
than once. 

To get such a basis one can take any spanning tree for the graph and then, for 
each non-tree edge, take the cycle formed by that edge and the unique tree-path 
connecting its endpoints. Finding a spanning tree for a graph is a classical computer 
science problem for which we use breadth- first search (see [8]). The spanning tree 
is also used for the task of finding a basis for M{P), as this requires finding paths 
in the Coxeter symbol from a fixed vertex to four others. Breadth-first search has 
the advantage of producing a short bushy tree which in turn gives short paths and 
cycles. 

7 Many computed examples 

In this section we present a number of examples including some unexpectedly com- 
mensurable pairs of groups. We also present some borderline cases of groups that 
are incommensurable, but have some of the invariants in common. 
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7.1 Lambert Cubes 



Recall from Section [5] that a Lambert cube is a compact polyhedron realizing the 
combinatorial type of a cube, with three disjoint non-coplanar edges chosen and 
assigned dihedral angles j, and ^ with l,m,n > 2, and the remaining edges 
assigned dihedral angles ^. Any reordering of {l,m,n) can be obtained by applying 
an appropriate (possibly orientation reversing) isometry, so when studying Lambert 
cubes it suffices to consider triples with I < m < n. 

In Table [1] we provide the invariant trace fields and the ramification data for the 
invariant quaternion algebras for Lambert cubes with small /, m, and n. 

It is interesting to notice that many of the above computations are done by 
hand in [13] , where the authors determine which of the "Borromean Orbifolds" are 
arithmetic by recognizing that they are 8-fold covers of appropriate Lambert cubes. 
(They call the Lambert cubes "pyritohedra." ) Our results are consistent with theirs. 

In this table and all that follow, we specify the invariant trace field kV by a 
canonical minimal polynomial p(z) for a primitive element of the field and, below this 
polynomial, a decimal approximation of the root that corresponds to this primitive 
element. 

We specify the real ramification data for the invariant quaternion algebra AT by 
a vector v whose length is the degree oi p{z). If the i-th root p{z) (with respect 
to Pari's internal root numbering scheme) is real, we place in the z-coordinate of v 
a 1 if AT is ramified over corresponding completion of kT, otherwise we place a 0. 
If the i-th root of p{z) is not real then we place a —1 in the i-coordinate of v. 

Below this vector v we provide an indication of whether AT ramifies at a prime 
ideal Vp over a rational prime p. In the case that AT ramifies at multiple prime 
ideals over the same p, we list multiple symbols Vp,T"p, . . ., etc. More detailed 
information about the generators of these prime ideals can be obtained in PARI's 
internal format using our scripts. 

The column labeled "arith?" indicates whether the group is arithmetic. In the 
case that the group satisfies conditions (1) and (2) from Theorems [13] and [14] but 
has elements with non-integral traces, we place a star next to the indication that the 
group is not arithmetic. (Some authors refer to such groups as "psuedo-arithmetic") 

7.2 Truncated cubes 

Another very simple family of hyperbolic refiection groups is obtained by truncating 
a single vertex of a cube, assigning dihedral angles j,^, and ^ to the edges entering 
the vertex that was truncated and ^ dihedral angles at all of the remaining edges. 
Since the three edges entering the vertex that was truncated form a prismatic 3- 
circuit, Andreev's Theorem provides the necessary and sufficient condition that 
J + ^ -|- i < 1 for the existence of such a polyhedron. In Table[2]we list the (/, m, n) 
truncated cubes for l,m,n < 6. 



20 



[I, m, n) 


kT 


disc 


yir ramification 


Arith'.' 


(3,3,3) 


X — X — X — X + 1 

-0.651387818 - 0.7587449567' 


-507 


[1,1,-1,-1] 



Yes 


(3,3,4) 


x' + l 

l.OOOOOOOOOOOOOOz 


-4 


[-1-1] 



No* 


(3,3,5) 


x^ -x'' - 3a;'' + - 3x^ ~x + l 
0.725191949 - 0.688546757i 


102378125 


[1.1,1,1,-1,-1,-1,-1] 



No 


(3,3,6) 


+ 5x^ - 3 
—2 35401 3863? 


-4107 


[LI, -1,-1] 

• ii • 3 


No* 


(3,4,4) 


- 2z^ - 2x + 1 
-0.3660254039 + 0.930604859i 


-1728 


[1,1,-1,-1] 




Yes 


(3,4,5) 


X** - x' - ISx" - ISx" + 95x* + 218x-' + 182a;^ + 7lx + 11 
-2.142201597+ 1.146040793i 


249761250000 


[1.1,1.1,-1.-1,-1.-1] 




No 


(3,4,6) 


x'^-x-' - Ux'^ + 33x - 6 
2.386000936- 1.441874268i 


-191844 


[1.1,-1.-1] 

• I' • 2 


No* 


(3,5,5) 


X** - 3x'' + x^ + 3x^ + - 31 + 1 
-0.2251919494 + 0.974314418i 


184280625 


[1,1,1,1,-1,-1,-1,-1] 



No 


(3,5,6) 


X** - x^ - 19x'' + 30x'^ + 741" - 72x^ - SlOx^ + 413x - 71 
—1.6937824749171233274 — 1.2086352565016507434? 


876653128125 


[1,1,1,1,-1.-1,-1,-1] 



No 


(3,6,6) 


X* - x^ - x^ - X + 1 
-0.651387818 - 0.758744956i 


-507 


[1.1,-1.-1] 



Yes 


(4,4,4) 


x* — x^ — 1 
0.786151377z 


-400 


[1,1,-1.-1] 



Yes 


(4,4,5) 


X** - 4x^ + 4x'' - 6x'^ + 19X-' - 14x'' + 4x^ - 6.T + 1 
-0.748606020 + 1.434441708i 


1548800000 


[1,1,1.1,-1.-1,-1,-1] 



No 


(4,4,6) 


x* - 2x'' - 6x + 9 
-0.822875655 + 1.524098309i 


-9408 


[1.1,-1.-1] 


No* 


(4,5,5) 


X** + bx*^ - 3x^ - 20x2 + 16 
1. 6754054322 


2714410000 


[1.1,1,1,-1,-1,-1,-1] 


No 


(4,5,6) 


x** - 3x'' - 25x'' + 95x" + SOx'-* - 3x-' - 1751x" + 2600x - 995 
-1.713666- 1. 9492541 


21059676450000 


[1.1,1.1,-1.-1.-1,-1] 




No 


(4,6,6) 


x* — X'' — x^ + 7x — 2 
1.280776406- 1.386060824i 


-10404 


[1.1,-1.-1] 

V2.V', 


No* 


(5,5,5) 


x" - x^ + x^ - X + 1 
0.809016994 - 0.587785252?: 


125 


[-1.-1,-1.-1] 




No 


(5,5,6) 


X** + 9x'' + Wx* - 27x2 + 9 
-2.027709945i 


5863730625 


[1,1,1,1,-1,-1,-1,-1] 


No 


(5,6,6) 


X** - 2x^ - 2x'' + 16x'^ - 15x'* - 19x'' + 43x2 - 22x + 1 
1.353181081 - 1.579017768i 


10637578125 


[1.1,1.1.-1.-1.-1,-1] 



No 


(6,6,6) 


x* — x~' — 3x2 _ ^ _^ I 
-0.895643923 - 0.4447718088i 


-1323 


[1.1,-1.-1] 



Yes 



Table 1: Commensurability invariants for Lambert Cubes. 



21 



{L m, n) 


kV 


disc 


AT ramification 


Arith? 


(3, 3, 4) 


x"^ — 2x^ + + 2x — 1 
1.207106781 - 0.978318343?: 


—448 


[1 1 -1 -11 
[J-, J-, J-, J-j 




Yes 


(3. 3, 5) 


- - 1 
0.786151377i 


-400 


fl 1 -1 -11 



Yes 


(3. 3. 6) 


+ 2x^ — 11 
2.112842071i 


—6336 


L-^i J-> J-J 



No* 


(3, 4i 4) 


+ 2 
1.414213562i 


_8 


[-1 -11 


Yes 


(3. 4, 5) 


T** + 1 4t'' + 57t'* -I- S6t^ + 41 
1.35712361i 


26869760000 


I— 1 —1 —1 —1 —1 —1 —1 —11 



No 


(3. 4. 6) 


x"* + 12x^ + 81 
1.224744871 - 2.738612788i 


57600 


[-1 -1 -1 -11 
[ J-, J-, J-j 




No 


(3, 5, 5) 


X* + 7x^ + 11 
-2.148961142?: 


4400 


[-1 -1 -1 -11 



No 


(3. 5, 6) 


_ 4^:7 _|_ 30^,6 _ g^j;-^' + 262x* — 384x'' + 978x^ — 684x + 1629 
1.3660254037844386468 - 1.7150676861906891827?: 


27,34871 040000 


[-1,-1.-1,-1.-1,-1.-1,-1] 
' ' ' 


No 


(3. 6, 6) 


x*^ — X + 4 

0.500000000 + 1.936491673i 


— 15 


[-1 -11 
1 -J, -i-j 


Yes 


(4. 4, 4) 


X-* - 2x'^ - 1 
0.643594253i 


— 1024 


n 1 -1 -11 
[1, 1, 1, ij 




Yes 


(4.4, 5) 


x* + 3x^ + 1 
1.618033989i 


400 


[-1 -1 -1 -11 



No 


(4,4,6) 


x-i + 2x^-2 
1.652891650i 


-4608 


[1,1.-1,-1] 



No* 


(4,5,5) 


X** - 2x*' - 8x^ - Sx" + 8x-' + 12x^ + 4.t - 1 
-0.707106781 + 0.544223826?: 


368640000 


[1,1,1,1,-1,-1.-1,-1] 


No 


(4,5,6) 


x» + 2x« - 39x'' - 130x2 - 95 
-1.05321208? 


-5042995200000 


[1,1.-1,-1.-1,-1.-1,-1] 



No 


(4,6,6) 


x-* + 2x^ + 4 

0.707106781 - 1.224744871? 


576 


[-1,-1.-1,-1] 



No 


(5,5,5) 


x-* - 5 

1.495348781i 


-2000 


[1,1,-1,-1] 



Yes 


(5,5,6) 


X** + Ox** - 13x'' - 66x2 _j_ 
2.466198614? 


12648960000 


[1,1.1,1.-1,-1.-1,-1] 



No 


(5,6,6) 


X-' + 21x2 ^ gg 

2.673181257? 


39600 


[-1,-1.-1,-1] 



No 


(6,6,6) 


x-* - 2x^ - 2x H- 1 
-0.36602.54039 + 0.930604859ii 


-1728 


[1,1.-1,-1] 



Yes 



Table 2: Commensurability invariants for some truncated cubes. 
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Figure 4: The Lobell polyhedron for n = 8. 



Figure 5: The "hexahedron" Hn- Edges labeled by an integer n are assigned dihedral 
angle ^ and unlabeled edges are assigned |. 

7.3 Lobell Polyhedra 

For each n > 5, there is a radially-symmetric combinatorial polyhedron having two 
n-sided faces and having 2n faces with 5 sides, which provides a natural general- 
ization of the dodecahedron. This combinatorial polyhedron is depicted in Figured] 
for n = 8. 

Andreev's Theorem provides the existence of a compact right-angled polyhedron 
Ln realizing this abstract polyhedron because it contains no prismatic 3-circuits or 
prismatic 4-circuits. 

An alternative construction of L„ is obtained by grouping 2n copies of the "hex- 
ahedron" shown in Figure [5] around the edge labeled n [30] . This construction is 
shown for Liq in Figure [9l We denote this polyhedron by H^, and note that, by 
construction, Ln and Hn are commensurable for each n. 

Of historical interest is that the first example of a closed hyperbolic manifold was 
constructed by Lobell [20] in 1931 by an appropriate gluing of 8 copies of Lq. (See 
also |30j for an exposition in English, and generalizations.) This gluing corresponds 
to constructing an index 8 subgroup of the reflection group in the faces of Lq , hence 
the Lobell manifold has the same commensurability invariants as those presented for 
Lq in the table above. It is also true that the Lobell manifold is arithmetic, because 
this underlying reflection group is arithmetic. Arithmeticity of the classical Lobell 
manifold was previously observed by Andrei Vesnin, but remained unpublished |28j . 

Furthermore, Vesnin observed in [29] that the if the Lobell polyhedron L„ is 
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kr 


disc 


AT ramification 


Arith? 




-x'-l 
0.7861511 


-400 


[1,1,-1,-1] 



Yes 


L, 


x' + 2 
1.414214 


-8 


[-1,-1] 


Yes 


Lj 


x'> + x'^-2x'^-l 
-0.6671H 


153664 


[1.1.-1.-1,-1,-1] 



No 


Ls 


x^ -2x^-1 
0.6435941 


-1024 


[1,1.-1,-1] 



Yes 


L9 


I** + 3x'' - 3 
1.5912541 


1259712 


[1,1,-1,-1,-1,-1] 



No 


Lio 


x^ + Sx^ + 1 
0.6180341 


400 


[-1,-1,-1,-1] 



No 


ill 


x"' + 4x** + 2x^ - 5x^ - 2.x2 + 1 
1.6378361 


-219503494144 


[1.1.1.1.-1,-1,-1.-1.-1,-1] 



No 




+ 2x^ - 2 
1.652892?: 


-4608 


[1,1.-1,-1] 



No* 


Lis 


+ 5a.io + 5,j.8 _ 5^.6 _ 7^.4 + 2x^ + 1 
1.664606?: 


564668382613504 


not computed 
(long computation) 


No 


Lu 


x'^ + X* - 2x^ - 1 
1.342363? 


153664 


[1,1,-1,-1,-1,-1] 


No 


Lu, 


X** + 5x* + 5x'' - 5x''^ - 5 
1.681396? 


-1620000000 


[1,1.-1.-1,-1,-1.-1.-1] 



No 


Lui 


X** + 4xS + 2x^ -4x^-1 
1.687530? 


-1073741824 


[1,1,-1.-1,-1,-1,-1.-1] 



No 


Ln 


not computed (long computation.) 








Lis 


x'' -3x^-1 
0.589319? 


419904 


[1.1.-1.-1,-1,-1] 


No 



Table 3: Commensurability invariants for the Lobell polyhedra. Note that L5 is the 
right-angled regular dodecahedron. 

arithmetic, then n = 5,6, 7, 8, 10, 12, or 18. He shows that the reflection group 
generated by L„ contains a (2,4, n) triangle group which must be arithmetic if -L„, is 
arithmetic, and apphes the classification of arithmetic triangle groups by Takeuchi 
|26j . In combination with our computations, Vesnin's observation yields: 

Theorem 15 The Lobell polyhedron Ln is arithmetic if and only if n = 5,6, or 8. 

Table [3] contains data for the first few Lobell polyhedra L„ and, consequently, 
the first few "hexahedra" Hn- 

It is interesting to notice that L7 and L14 have isomorphic invariant trace fields, 
which are actually not the same field (one can check that the specified roots generate 
different fields). This alone suffices to show that L7 and L14 are incommensurable. 
Further, albeit unnecessary, justification is provided by the fact that their invariant 
quaternion algebras are not isomorphic, since AT{Lj) has no finite ramification, 
whereas AT{Li/C} is ramified at the two finite prime ideals Vt,V'-j lying over the 
rational prime 7. 
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q 


kV 


disc 


Arith? 


7 


;r" - 2.x'' + a--* - Ax-' + 3x'^ + 3.x - 1 
-0.400969 - 1.444370? 


-199283 


Yes 


8 


- 4x^ + 6x''' - 7x* + 2x^ + 4x - 1 
1.423880 + 1.4948381 


473956352 


No 


9 


_ 4^.9 ^ 27x» - 6a:^ - 26x'' + Gx'^ - 39.x* + 34i-» - ISx^ + 6x - 1 
-1.646799- 1.864938i 


5879193047138304 


No 


10 


- 4x^ + 7x« - 2x^ - Ss" + 2x^ + 2x-2 + 4x - 4 
1.451057+ 1.5538934 


-380000000 


No 


11 


x"' - 3x» + x^ + x^ + 9x« - 18x^5 + 5x* + 9x-* - 6x^ - x + 1 
-1.114354- 1.200301Z 


-14362045027 


No 


12 


X* - 4x^ + 6x'' - 4x^ - 9x'' + 20x'' - lOx^ + 1 
0.500000 + 1.7018414 


514916352 


No 


13 


xi2 - 5x" + lOx^" - 17x8 + 32x** - 34x' + 15xS - llx^ + 2x'' + 22x'^ - 9x^ - 6x + 1 
-0.470942 + 1.5965981 


-214921388792591 


No 


14 


x^'^ - 62;" + 16x"' - 18x" - 2.x** + 22x' - 8x'' - 8x"^ - 16.x* + 18x'' + 9x^ - 8x + 1 
1.474928- 1.6051114 


-581980365811712 


No 


15 


x» - 2x7 _ 2^1.6 _|_ 53.5 _ 5^.4 _ 4^.3 ^ 3^2 ^ 3,,. ^ I 
-0.309017 - 0.336995« 


183515625 


No 



Table 4: Commensurability invariants for truncated prisms. 

7.4 Modifying the Lobell 6 polyhedron 

Let P be a compact hyperbolic polyhedron with all right dihedral angles and at least 
one face F with 6 or more edges. If one (combinatorially) splits the face F into two 
faces Fi and F2 along a new edge e, the resulting polyhedron can be realized with 
all right dihedral angles so long as both of the faces Fi and F2 have 5 or more edges. 
The inverse of this procedure is described in |17] . 

In Figure[6]we illustrate various right-angled polyhedra obtained from the Lobell 
6 polyhedron by adding such edges and their commensurability invariants. 

7.5 Truncated prisms 

Maclachlan and Reid consider triangular prisms with dihedral angles as labeled on 
the left hand side of Figure \7\ Aside from the tetrahedral reflection groups and 
some of the Lambert cubes, this family of prisms is one of the few cases that can be 
computed "by hand." Natural candidates for testing our program are the truncated 
versions shown on the right hand side of Figure [71 

In Table U] we show the arithmetic invariants that we have computed for the 
truncated prisms for q = 7, ... ,15. Each of the quaternion algebras is ramified at 
all real places and not at any finite places, so we omit ramification data. 

7.6 Doubly-truncated prisms 

For q = A and 5, there exist compact polyhedra realizing two doubly-truncated 
prisms pictured in Figure El 
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kr -.x^ -x^ + 3x* + 4a;2 + 1, 0.762867j kV : x^ - 3x'^ + 5x^ - 7x^ + 6x* - Qz^ + 5x^ - 4x + 1 

disc = 36096064 

Ar : [-1,-1,-1,-1,-1, -1,-1,-1], 



disc = 36096064 0.167273 - 0.768150i 

Not arithmetic '^^^'^ = "62410000 



: [1,1,-1,-1,-1,-1,-1,-1],0 
Not arithmetic 



Figure 6: Right angled polyhedra obtained by adding edges to the Lobell 6 polyhedron. 
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Figure 7: On the left is the family of prisms considered in [22] . On the right is the family 
considered in this section. As before, edges labeled by an integer n are assigned dihedral 
angle ^ and unlabeled edges are assigned dihedral angle |. 




Figure 8: Edges labeled by an integer n are assigned dihedral angle ^ and unlabeled edges 
are assigned dihedral angle |. (In particular, the triangular faces are at right angles to 
each adjacent face.) 
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When g = 4, the invariant trace fields are equal, both equal to Q(a), with a 
is an imaginary fourth root of 2, and ATi = AT 2 — ^ q(q) ^ ^ • Indeed, one can 
show this using the fact that Fi and r2 contain subgroups isomorphic to ^4, see the 
comment at the end of Section [3l However, Fi and r2 are incommensurable since 
Fi is arithmetic, while F2 has non-integral traces. 

If we repeat the calculation with q = 5 both groups have invariant trace field gen- 
erated by an imaginary fourth root of 20 and have isomorphic invariant quaternion 
algebras. However, both groups have non-integral traces, so neither is arithmetic. 
Thus, we cannot determine whether these two groups are commensurable, or not. 
(See Subsection 17. 8[ below, for similar examples of non-arithmetic polyhedra with 
matching pairs (kT,AT).) 

There is a good reason why these pairs have the same invariant trace field: 
The invariant trace field for an amalgamated product of two Kleinian groups is the 
compositum of the corresponding invariant trace fields. See Theorem 5.6.1 from |22j . 
The group on the left can be expressed as an amalgamated product obtained by 
gluing the "top half," a singly truncated prism to the "bottom half," another singly 
truncated prism (that is congruent to the top half) along a {q, q, g)-triangle group. 
The group on the right can be expressed in the same way, just with a different gluing 
along the {q, q, q) triangle group. This construction is similar to the construction of 
"mutant knots," for which commensurability questions are also delicate. See page 
190 of [22]. 



7.7 Unexpected commensurable pairs 

As noted in Section [3l one way to find unexpected pairs of commensurable groups 
is to verify that two groups are arithmetic, have the same invariant trace field, 
and have isomorphic invariant quaternion algebras. In this section, we describe the 
commensurability classes of arithmetic reflection groups in which we have found 
more than one group. It is very interesting to notice that Ian Agol has proven that 
there are a finite number of commensurability classes of arithmetic reflection groups 
in dimension 3, [2], however there is no explicit bound. 

The reader may want to refer to the Arithmetic Zoo section of [22] to find other 
Kleinian groups within the same commensurability classes. 

It is easy to see that the dodecahedral reflection group is an index 120 and 60 
subgroup of the reflection groups in the tetrahedra T2 and respectively. However, 
for most of these commensurable pairs the commensurability is difficult to "see" 
directly by finding a bigger polyhedron Q that is tiled by reflections of each of the 
polyhedra within the commensurability class. 

For example, the commensurability of the right-angled dodecohedron L5 with 
the (4, 4, 4)-Lambert cube was first discovered using our computations and only 
later did the authors find the explicit tiling of the dodecahedron by 8 copies of the 
(4, 4, 4)-Lambert cube shown in Figured! We leave it as a challenge to the reader to 
observe the commensurability of the (3, 3, 5)-truncated cube with the dodecahedron 
in a similar way. (Note that at least 13 copies of the dodecahedron are required!) 
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kT 


disc 


AT finite ramification 


+ 2, root: 1.414213562z 
The (3,4,4) truncated cube 
The Lobell 6 polyhedron 


-8 


vol ^ 1.0038410 
vol ^ 6 ■ 1.0038410 


- x2 - 1, root: 0.78615Z 
Compact tetrahedron from p. 416 of [22] 
Compact tetrahedron T4 from p. 416 of [22] 
The (4,4,4) Lambert cube 
The (3,3,5) truncated cube 
The Lobell 5 polyhedron (dodecohedron) 


-400 




vol ^ 0.03588 
vol ^ 2 ■ 0.03588 
vol ^ 15 ■ 0.03588 
vol ^ 26 ■ 0.03588 
^ 120 ■ 0.03588 


X* - 2x-^ + x^ + 2x-l, root 1.207106781 - 0.978318343i 
The prism on the left of Figure [7] with g = 8 
The (3,3,4) truncated cube 


-448 




vol ^ 0.214425456 
vol 16 • 0.214425456 


- - - X + 1, root: -0.651387818 - 0.758744956z 
The (3,3,3) Lambert cube 
The (3,6,6) Lambert cube 


-507 




vol ^ 0.324423449 
vol fa f ■ 0.324423449 


- 2x2 - 1, root: 0.643594253? 
The (4,4,4) truncated cube 
The Lobell 8 polyhedron 


-1024 




vol ^ 1.1273816 
vol ^ 8 ■ 1.1273816 


x^ - 2x3 _ 2x + 1, root: -0.366025 + 0.930605z 
The (3,4,4) Lambert cube 
The (6,6,6) truncated cube 


-1728 




vol ^ 0.4506583058 
vol ^ 6 ■ 0.4506583058 



Table 5: Commensurablity classes of arithmetic reflection groups. 
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Figure 9: Tiling the Lobell 10 polyhedron with 10 copies of the (4,4,5) truncated cube. 



For n = 3 and 4 we see the commensurability between the Lobeh 2n polyhedron 
and the (4, 4, n) truncated cube. This commensurabihty is a general fact for any 
n: one can group 2n of the (4, 4, n) truncated cubes around the edge with label 
n forming the right-angled 2n Lobell polyhedron. This construction is shown for 
n = 5 in Figure [9l 

The approximate volumes were computed using Damian Heard's program Orb 
(Orb provides twice the volume of the polyhedron, which is the volume of the 
orbifold obtained by gluing the polyhedron to it's mirror image along its boundary.) 

7.8 Pairs not distinguished by kT and AT 

We found three pairs of non-arithmetic polyhedral reflection groups which are indis- 
tinguishable by the invariant trace field and the invariant quaternion algebra. The 
first of these is the pair of doubly truncated prisms with q = 5, from Subsection 17.61 
Each of these has the same volume, approximately 2.73952694 since the second can 
be obtained from the first by cutting in a horizontal plane and applying a | twist. 
We leave it as an open question to the reader whether the two polyhedra in the first 
pair are commensurable. 

The second and third pairs are the (4, 4, 5) truncated cube and the Lobell 10 
polyhedron and the (4, 4, 6) truncated cube and the Lobell 12 polyhedron, each of 
which fit into the general commensurability between the (4, 4, n) truncated cube 
and the Lobell 2n polyhedron, as described in the previous subsection. 
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A Computing finite ramification of quaternion 
algebras 

Recall from Section S] that to determine whether two quaternion algebras over a 
number field F are isomorphic it suffices to compare their ramification over all real 
and finite places of F. 

To compute the finite ramification of a quaternion algebra A = (j^^ we first 
recall that there are a finite number of candidate primes over which A can ramify: 
A necessarily splits (is unramified) over any prime ideal not dividing the ideal {2ab) . 
To check those primes V dividing {2ab), we apply Theorem |4] which states that A 
splits over V if and only if there is a solution to the Hilbert Equation aX'^ + bY'^ = 1 
in the completion Fiy. Here u is the valuation given by u{x) = c"^(^') on Rf, see 
Section HI 

Because it is difficult to do computer calculations in the completion F^, where 
elements are described by infinite sequences of elements of F, we ultimately want 
to reduce our calculations to be entirely within Rp. Hensel's Lemma, see |18] . is 
the standard machinery for this reduction. Happily, for the problem at hand, the 
proper use of Hensel's Lemma has previously been worked out, see [9] (whose authors 
write ly-p instead of our n-p and | • | instead of our i/). We use the same techniques 
they do but without many of the optimizations, which we find unnecessary for our 
program (this is not the bottleneck in our code). The following theorem provides 
the necessary reduction: 

Theorem 16 Let V be a prime ideal in Rf, let a,b € Rp be such that n-p {a), n-p [b) € 
{0, 1}, and define an integer m as follows: ifV\2, m = 2np{2) + 3 and ifV /2 then 
m = 1 ifnp{a) = n-p{b) = and m = 3 otherwise. 

Let S be a finite set of representatives for the ring Rp/V"^. The Hilbert Equation 

aX'^ + bY'^ = 1 

has a solution with X and Y £ F^ if and only if there exist elements X\Y\ and 
Z' £ S such that 

viaX''^ + _ ^'2^ < 
and max{i/(X'),i^(y')>i^(^')} = 1- 

Recall that < c < 1 is the arbitrary constant that appears in the definition of v. 
Also note that the condition that np{a),np{b) G {0, 1} is no real restriction since 
one can always divide the elements of F appearing in the Hilbert symbol by any 
squares in the field without changing A. 

Our Theorem[T6]is a minor extension of Proposition 4.9 of [9]. Their proposition 
only applies to dyadic primes, corresponding to m = 2n-p(2) + 3. The other two 
cases can be proved analogously using, instead of their Lemma 4.8, the following 
statement: 
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Lemma 17 Suppose that v is the valuation corresponding to a non-dyadic prime 
V. Let X,X' be in Fy and suppose that v{X) < 1 and v{X — X') < for some 
non-negative integer k. Then, v{X'^ — X'"^) < c?^ . 

The reason we make this minor extension of Proposition 4.9 from [9] is that we 
do not implement the optimizations for non-dyadic primes that appear in SNAP, 
instead choosing to use the Hilbert Equation in all cases. 

Thus, the problem of determining finite ramification of A reduces to finding solu- 
tions va. Rp to the Hilbert inequality from Theorem [TBI Our program closely follows 
SNAP [To], solving the equation by exploring, in depth-first order [8], a tree whose 
vertices are quadruples (X, y, Z, n), where v{aX'^ -\- hY"^ — Z'^) < c". The children 
of the vertex {X,Y, Z,n) are the vertices {Xi,Yi, Zi,n + 1) where Xi = X,Yi = 
Y, and Zi = Z (mod P"). In fact, the condition max{z^(X), z^(y), z^(Z)} = 1 
means that the search space is considerably reduced because there must be a solu- 
tion to the inequality with one of X, Y, Z equal to 1. Indeed, if I'iY) = 1 then Y is 
invertible modulo any power of V and {XY~^, 1, ZY~^) is also a solution. So, our 
program only searches the three sub-trees given by fixing X,Y, or Z at 1. In the 
case that we search all three trees up to level n = m unsuccessfully, then Theorem 
[TBI guarantees that there is no solution to the Hilbert Equation, and consequently 
A ramifies over V. 
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